1 



Extremal dependence: some contributions 



Helena Ferreira Department of Mathematics, University of Beira Interior, Covilha, Portugal 

Marta Ferreira Department of Mathematics, University of Minho, Braga, Portugal 

Abstract: Due to globalization and relaxed market regulation, we have assisted to an increasing of ex- 
tremal dependence in international markets. As a consequence, several measures of tail dependence have 
been stated in literature in recent years, based on multivariate extreme-value theory. In this paper we 
present a tail dependence function and an extremal coefficient of dependence between two random vectors 
that extend existing ones. We shall see that in weakening the usual required dependence allows to assess 
' the amount of dependence in d-variate random vectors based on bidimensional techniques. Very simple 

estimators will be stated and can be applied to the well-known stable tail dependence function. Asymptotic 
bJQ' normality and strong consistency will be derived too. An application to financial markets will be presented 

^ ' at the end. 
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^ ' 1 Introduction 

I Dependence between extremal events have increased in recent time periods in financial markets, 

especially during bear markets and market crashes. The globalization and the lack of supervision 
are well-known contributions for this phenomena. Therefore, modern risk management is highly 
interested in assessing the amount of extremal dependence. The concept of tail dependence is the 
current tool used to this end, although it was first introduced far back in the sixties (Sibuya |25) . 
1960; Tiago de Oliveira 1962/63). Tail dependence coefficients measure the probability of 
' occurring extreme values for one random variable (r.v.) given that another assumes an extreme 

^\ . value too. These coefficients can be defined via copulas of random vectors which refers to their 

dependence structure concerning extreme events independently of their marginal distributions. 
. The tail dependence coefficient, 

o . 

A = lim P{Fx{X) > 1 - t\FY{Y) > 1 - i), (1) 

where Fx and Fy are the distribution functions (d.f.'s) of X and Y, respectively, is perhaps the 
, most referred in literature and characterizes the dependence in the tail of a random pair (A, Y) , 

^ ' i.e., A > corresponds to tail dependence and A = means tail independence. There are several 

^ references on this topic (besides the two above) and thus we point out only some of them: Ledford 

and Tawn ([HI US], 1996, 1997), Joe ([TO], 1997), Coles et al. ([2], 1999), Embrechts et al. (0, 
2003). 

Multivariate formulations for tail dependence coefficients can be used to describe the amount 
of dependence in the orthant tail of a multivariate distribution (Marshall-Olkin [17] 1967; Wolff 
[H] 1980; Nelsen [H] 1996; Schmid and Schmidt jH] 2007; Li [ll[T3[TO] 2006, 2008, 2009, among 
others). These have been increasingly used in the most recent and higher demanding times. Most 
of the multivariate measures consider that extremal events must occur to all the components of the 
random vector, and obviously they are more complicated to deal with and to understand than in 
the bivariate case. Not surprisingly, applications hardly go any further than the three-dimensional 
case. 

But maybe this is a too demanding condition and the occurrence of at least one extremal event 
in sub-vectors of a random vector can be enough to assess dependence. As already mentioned. 
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financial markets are increasingly connected and the occurrence of at least one market crash, for 
instance, in Europe, will certainly influence a negative behavior in USA markets. 

Based on this, we define a new tail dependence function for a random vector as a measure of 
the probability of occurring extreme values for the maximum of one sub- vector given that the max- 
imum of another assumes an extreme value too. At the unit point, this function gives rise to the 
here called extremal coefficient of dependence as it relates to the well-known extremal coefficient 
(Tiago de Oliveira 1962-63, Smith 1990). These extend, respectively, the concept of upper tail 
dependence function and upper tail dependence coefficient already stated in literature (see Schmidt 
and Stadtmiiller (2006) and references therein). In deriving the moments of the random variables 
involved in this approach, we find very simple estimators that can be also applied to the well-known 
stable tail dependence function (for a survey on this function see e.g. Beirlant et al. [T], 2004). 
Asymptotic normality and strong consistency are proved. 

This paper is organized as follows. In Section[2]we define our new upper-tail dependence function 
and the extremal coefficient of dependence. We present some properties and examples. We also 
analyze the case of asymptotic independence. In Section |2] we present estimators and derive the 
respective properties of asymptotic normality and strong consistency. Section 2] illustrates our 
approach through an application to financial data. 



2 Extremal dependence between two random vectors 

Let X = {Xi, Xd) be a random vector with d.f. F and continuous marginal d.f.'s Fi. For 
/ C {!,..., d}, define M(J) = \J Fi{Xi) and X/ the sub-vector of X having r.v.'s with indexes 
in /. Consider Cp the copula function of F, i.e., 

F{xi,...,Xd)=CF{Fi{xi),...,Fd{xd)). (xi,...,Xd) e M'^. (2) 

We are going to study the dependence between extremal events concerning two sub-vectors, 
X/j and X/2, where Ii and I2 are disjoint subsets of {1, d}. 

We start by extending in Definition 12.21 the concept of upper tail dependence function (see 
Schmidt and Stadtmiiller (2006) and references therein) and from this we define a new tail depen- 
dence coefficient between two random vectors. 

Definition 2.1. Let Ii and I2 he two non-empty subsets 0/ {1, c?}. The upper-tail dependence 
function of~X.j-^ given X/^ is defined as, for {x,y) G {0,oo)'^, 

M(/2)>1-|), (3) 

provided the limit exists. 

By taking x = y = 1, we have 

Mih)>l-j), (4) 
which is a tail dependence coefficient greater than the one considered in Li and Sun (2008), 

7 = /nn P( fl F,{Xi} > 1 - ^1 U ^'(^^) > 1 - i)' 
ieii ieh 

which in turn is greater than the coefficient of Li (2009) for Ii — {1, d} — I2, 



A\j'^''\x,y) = lim p(M{h) > I-7 



Aij'^^'Hl, 1) - lim p(m{Ii) > 1 - 

t— !-00 \ t 



The tail dependence coefficient a[|^'^^^(1, 1) give us information about the probabihty of occur- 
ring some extreme value in {Fi{Xi)^i G /i} given that some extreme value occurs in {Fi(Xi)^i € 

Before presenting the properties of function A[|^ (x, y) that will be the basis for the definition 
of our coefficient, consider the following notation: 

for {x,y) € (0,oo)2, C Q {1, and i € {1, ...,d}, let 

a\'"''\x,y)=xli,{t) + yli,{t) + ^lT^it), (7) 

where l(-) is the indicator function, and 

l^'^''^\x,y) = -\ogF{a[''''^\x,y), ...,a[,''^'^-\x,y)), (8) 

with the convention that, when some of the arguments of F are 00 we understand the limit of F 
as those arguments tend to 00. 

If is a multivariate extreme value distribution (MEV) with unit Frechet marginals, we have 

^(/ij2)(2.-i^.^-i) = -log(cxp(-a:;))^^i-^2 = xei.^i,, (9) 
where e/iU/2 the extremal coefficient of X/ju/2 (Tiago de Oliveira 1962-63, Smith 1990). 

Proposition 2.1. If F is a MEV distribution with unit Frechet marginals, then function A^^^^'^\x, y) 
is defined and verifies 

A^J^^'^\x,y) = l+^-^-'''''''^^''^y~'\ 
Dem. We have 

A(^^I^-)U- ,) hn. 1 I ^ - ^(^^('^^ ^ ^ - ^) ^ ^ ' - ^ ^ - ^) (11) 

^^^V)- }^^+ ^_p^M{h)<l-^) 1 - P(M(/2) < 1 - f ) • ^ ^ 

On the other hand, for a'^^'^^\u,v) = ulj^{i) + vli2{i) + ljjj]7j(i), it holds 
^fim - t logP(A/(/i) < 1 - f , M(/2) < 1 - f ) 

= lijn - logC.(a(^-^^'((l - f )*, (1 - f )*),..., 4^-^^) ((1 - f (1 - f 

= - log Cf (a['''^^^ ( exp(-x), exp(-y)) , a^^'"^'-* ( exp(-x), exp(-2/))^ 

= - logF(aJ^-^^' {x-\y-') , a^J'^'^'> {x'^y-'] 
- l('^^'-\x~\y~^). 

Therefore, dividing the numerator and denominator of the fractions in pip by t, we obtain 



(10) 



(12) 



;(0J.)(y-l,y-l) ;(0,/.)(y-l^y-l) 

- log(exp(-a;))'^^i l^^^'^^^x'^ ,y^'^) 



(13) 



log(cxp(-y))'^2 - log(exp(-?/))^^2 



□ 
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Therefore, under the conditions of Proposition 12. 1[ we have 
and we wih denote this common value as A^l^'^^-* (x, y). 

Definition 2.2. The upper-tail dependence function for random vector (X/^jX/^) with d.f. MEV 
and unit Frechet marginals is defined as 

A[j'^'-^\x,y)^xej,+yej,-l^'^-'^\x-\y-') (14) 

and the extremal coefficient of dependence between X/^ and X/^ is given by A\j''^^\l,l), which we 
denote e(/i,/2) ^^'^ hence 

^{iiM = + - e/iU/2- (15) 

The upper-tail dependence function p4l) generalizes the relation of Huang (1992) corresponding 
to h = {1} and h = {2}, 

kuix.y) ^ X + y - lfp^^x^),F^(x^)){x^y)^ (16) 
where the stable tail dependence function in the right-side is given by 

l(F,ix,),F,ix,))i^.y) = > 1 - 7 V F^{X^) > 1 - f ). (17) 

Observe that by pTjl we also obtain 

i<^'-''-\x-\y-^) - \\mtp(M{h) > 1 - ^ V M{h) > 1 - f ) = ifM(i,)Mim{^^y)- (18) 

Moreover, the upper-tail dependence function in (jl4p can be can be viewed as an extension of the 
bivariate upper-tail dependence function of Schmidt and Stadtmiiller ([H], 2006), defined as 

AfF,(X,).i=^.(X.))(^,y) = /ij^^^(^l(^l) > 1 - f >^2(X2) > 1 - f ), (19) 

by taking in this limit the random pair {M{Ii),M{l2)) instead of {F{Xi),F{X2)). At the unit 
vector, the Schmidt and Stadtmiiller upper-tail dependence function corresponds to the tail de- 
pendence coefficient A in ([1]), i.e., A = ^fpiiXi) F2(X2))^^^ 

In the following, we present the expression of the tail-dependence function A[|^'^^^(a;, y) and 
the value of the corresponding extremal coefficient ^^i-^.i^) for a d-variate random vector X with 
well-known distribution functions for its margins. 



Example 2.1. Consider vector X with unit Frechet margins and copula function Cx(ui, Ud) = 
nl^i rife^-oo "^V"^ ^ ■■■ ^ w^""*' where Uj £ [0, 1], j = 1, d, and {azfcj, — oo < k < oo, 1 < j < 
d,l > 1} is a family of non negative constants such that X^i^i X]fc°=-oo '^I'^j ^ 1> i = li d,- The 
distribution of X is the MEV marginal distribution of multivariate maxima of moving maxima 
processes considered in Smith and Weissman (\24^, 1996). We have 



oo oo 



/(^-^=)(x,y) = -logC(e-'^i"(--^),...,e-'^^"(-^^)) ^ y a-\x,y)aik,. 

l — \ k— — oo j — 1 

Therefore, 

A[j^^'^\x,y)=l^^^''\x-\x~')+l^'^'^\y-\y-')-l^'^'^^\x-\y-') 

OOOO OOOO OOOOy \ 

1=1 k=-cx} jeh 1=1 k=-oc jel2 1=1 k=-oa ^ j6/i 
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and 

OO CXD OO CXD CXD CXD 

/ — 1 — OO je/i l — lk— — ooj^l2 / — 1 fc— — OO je/iU/2 

Illustrating with 

/i = {1, 2} and I2 ~ {3, 4}, we obtain 

^\J"''\^.y) = (I + 1 + 1 + i)- + (I + 1 + l)y - ((4 V yl) + (4 v 4) + (4) + (4 v 4)) 



-y - ((^1 V yi) + (xf V 4) + xl + (xi V 



9 I 14 / 1 I Z I 2 I 6\ _ 7 

v8 ' 8 ' 8 ^ 8/ 8 



_ 9 I 14 (i 

Similarly, if Ii = {1,2} and I2 = {4} we obtain 



and 



k['y^\x, y)^lx + y- ({x\ V 4) + (4 V 4) + x| + {x\ V yf)) 



e(/i,/2) - 8+ ^^ 18 + 8 + 8 + ij- 8- 



Example 2.2. For i/ie symmetric logistic model we have 

d ^ 

mt/i G (0, 1], x,y > 0. Therefore, 

A^J''^'\x,y) = + 

j&ii j&h jell jeh 

= \hf^ + \l2\'y^{\h\x'/' + \l2\y'n' 

and 

eiii.i2) = \hf + \l2\'-{\h\ + \l2\f. 



Proposition 2.2. Under the conditions of Provosition we have 
ft) 0<A^jj"''\x,y)<xei,Ayei, 



(ii) < e(/i,/2) < '^h ^ 
Dem. 



(i) The left inequality is straightforward by the definition of A^^^'^^' (cc, y) in (fH)l . Observe also 
that, since X has MEV distribution, it is associated (in the sense of Joe [10], 1997; Theorem 
6.7) and hence, for all /i,/2 C {l,...,d}, 

F{ay'^\x-\y~^),...,a'-y''\x-\y~^)) 
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leading to the same conclusion, i.e., 
On the other hand, 

< F{a['^''\x-\x-^),...,a^J^^'\x-\x-'))AFia^^''^\y-\y-^),...,^ 
and hence 

< ?(^i-0)(.T-i,x-i) + /(0>^^)(y-i,y-i)- (/(^i*(x-i,.T-i) V/(0-^^)(y-i,j/-i)) 
= xe/i A ye/2 . □ 

The result in (i) agrees with the one for the bivariate case. Observe that, from the proof above 
we can also conclude that the boundary cases correspond to, respectively, independence and total 
dependence. 

Remark 2.1. With the conventions 1/0 := oo and l/oo ;= 0, we can define A\j^''^^\x,y) in 
[0,oo]2\{(oo,oo)} and found A\j''^'\0,y) = = A\j''^'\x,0), A\j''^'\oo,y) = yei, and A\j''^'\x,oo) = 

Proposition 2.3. Under the conditions of Provosition [K7\ and Remark \2.1[ for each y > 0, the 
partial derivative dA^^^'^^^/dx exists for almost all x > 0, and 

Q<i^A\^'''\x,y)<\h\. 

Similarly, for each x > 0, the partial derivative dA^^^'^^^/dx exists for almost all y > 0, and 

0<^A'u'-''\x,y)<\l2\. 

Also, the functions x i— ^ dA^j^'^^"^ (x, y)/dy and y dA^'^^'^^^x, y)/dx are defined and non decreas- 
ing almost everywhere on [0,oo). 

Dem. The function A[|^'^^' (x, y) is 2-increasing since a bivariate d.f. is 2-increasing. By 
Remark |2. II we conclude that A^^'^^\x,y) is grounded. Hence, applying Lemma 2.1.5. in Nelsen 
(dg, 2006) we have, for {x,y),{x* ,y*) € [0, cx)]2\{(oo, oo)}, 

|A[/-^^)(x-,y)-A[/-^^'(x*,y*)| 

< \mit(\P{M{h) > 1 - f ) - P(M(/i) > 1 - ^)| + \P{M{h) > 1 - f ) - P{M{h) > 1 - V)l) 

< |/i||x-a;*| + |/2||2/-y*|. 

Now, the proof is straightforward from Theorem 3 in Schmidt and Stadtmiiller ([Hj, 2006). □ 

Remark 12.11 and Propositions 12.21 fi) and 12.31 extend, respectively. Theorems l.i), 2.i) and 3 
of Schmidt and Stadtmiiller ([H], 2006). Moreover, given the above mentioned relation between 
A\j^'^^\x,y) and the bivariate upper-tail dependence function A^p.^^.^^^ p^^-^^-^-^{x,y) in (fTO)) . the 
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properties ii)-v) of Theorems 1 and 2 of Schmidt and StadtmiiUer (|22j, 2006) are straightforward 
for A[/-^^)(x,y). 



We now discuss the case of tail independence between M{Ii) and M{l2) and hence extend our 
context beyond a MEV distribution. 

Notice that, in case of tail dependence between r.v.'s Fi{Xi) and F2{X2), the mapping 

iK^p(^^i(Xi)>l-|,F2(X2)>l-|) (20) 

is regularly varying of order —1 at oo, and so an homogeneity property holds for large t. However, 
if (Fi (Xi), ^2(^2)) is tail independent, this latter does not hold and an adjusted homogeneity 
property can be obtained by assuming that (|20p is regularly varying of order — l/ryatoo, ?7<1 
(the case rj = 1 corresponds to tail dependence). Coefficient 77 is the coefficient of tail dependence 
introduced in Ledford and Tawn (1996, 1997). 

Thus being, if we assume that (|20p is regularly varying of order — 1/r/ at 00, i.e., 

P{F^{X,)>l~x/t,F2{X2)>l-y/t) 

hm -, ; ; ; ; ^ — = C (X , v) (21) 

t^<x,p(Fi(Xl)>l-lA,F2(X2)>l-l/t) ^ '^^ ^ ' 

for {x, y) £ [0, oo)^, where c* is homogeneous of order I/77 for some 77 € (0, 1] and c*(l, 1) = 1, then 
t 1-^ P(^Fi{Xi) > 1 — 1/^,^2(^2) > 1 — 1/t) is regularly varying at 00 with index —l/rj (choose 
X = y in (|2ip ). and hence we can write 

P(Fi(Xi) > 1 - l/t,F2{X2) > 1 - l/t) = t-^/'''L{t) (22) 

where L is a slowly varying function at 00 (i.e., L{tx)/L{t) — >■ 1, as t — >■ 00, for any x > 0). 
Observe that rj dominates the speed of convergence of P(Fi(Xi) > 1 - 1/^,^2(^2) > 1 - l/t) to 
0. If 77 < 1 then Fi{Xi) and ^2(^2) (and thus Xi and X2) are asymptotically independent (or tail 
independent). In this case, the tail dependence coefficient A in ^ is null. Conversely, asymptotic 
dependence holds if 77 = 1 and L{t) a > 0, as i — > 00, and we have A > 0. If 77 = 1/2 we 
have (almost) independence (perfect independence if L{t) = 1 and ([?!]) holds with c*{x,y) = xy). 
The cases 77 S (0, 1/2) and 7y € (1/2, 1) correspond to asymptotically negative independence and 
to asymptotically positive independence, respectively. Roughly speaking, coefficient governs a 
kind of a pre-asymptotic tail behavior that allows to better estimate the probability of extreme 
events in case of tail independence. A bivariate extreme value distribution (BEV) allows only tail 
dependence (77 = 1) or independence (r/ = 1/2), since 

P(Fi(Xi) > 1 - llt,F2{X2) > 1 - l/t) ^ (2 - /»i>'™(l, l))/t + ((/(ti>'™(l, 1))V2 - l)/t^ 

as t -> 00. For a discussion on this topic see, for instance, Ledford and Tawn ([T^, 1996), Draisma 
et al. (O, 2004) and Drees and Miiller (g], 2008). 

Now assume that (PT|) holds for random pair (A/(/i), A/(/2)), i.e., 

P{M{h)>l-x/t,M{l2)>l~y/t) 
}^2o P{M{I,)>l-l/t,M{l2)>l-l/t) = 

for {x,y) £ [0,03)^, where c^j-^ j,) is homogeneous of order l/r]{^i^j^) for some ViiiJ^) ^ (0' 1] ^^^^ 
C(/i,/2)(l. 1) = 1- Taking x ^y in (UHl), one obtains that P{M{Ii) >l-l/t, Ml^h) > 1 - l/t) is 
regularly varying at 00, i.e., 

P(M(/i) > 1 - l/t.Mih) > 1 - l/t) = t-i/''(^i>^^)P(,^,,,)(t), (24) 
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where L(^i^j^'^{t) is a slowly varying function at oo. Coefficient is now a measure of the 

speed of convergence of P{M{Ii) > 1 — A/(/2) > 1 — to and is, therefore, a coefficient 
of tail dependence between M(/i) and M{l2), with analogous conclusions derived for rj above. 
Similarly, in a MEV we obtain, as < — cx), 

P{M{h > 1 - 1/t, M{h) > 1 - l/t) ^ {ei, + ej, - ei,ui.)/t + {el^j^ ~ e\ - eX)/{2t^). 

Hence it only occurs asymptotic dependence whenever ^[i-^j^) ~ e/i + e/2 ^ ^huh > (with 
^(^1,^2) = 1)1 otherwise independence (?7(/i,/2) = 1/2). 

In the next result we compute ?7(/i./2) and found that it is given by the maximum coefficient 
Vi G h,j € h- 

Proposition 2.4. Suppose that {24^ holds and 

P( min (F,iX,),F,{X,)) > 1 - l/t) = i-i/^-'L,,, , (t) (25) 
holds for all ^ ^ I <Z Ii and 7^ ,/ C /2, where L^u , is a slowly varying function at 00. Then 

Dem. First observe that if /' C / and J' C J then 

1 > t-i/'^.',.'L,^, ^,(t) > t-'^^^'^-'L,,Jt). (26) 

We have that 

P(V,e,, F.,{X.,) > 1 - l/t, V,e/2 P^X,) > 1 - l/t) 
= Pi{J^^I^F.^{X^) > 1 - l/W,eiAP^(^j) > 1 - l/i}}) 

= E0^sc/,(-l)'^'+'^(ae5{^»(^O > 1 - l/i},Ue/2{^.;(^^) > ^ " VO) (27) 

— V V MS'l+|T|+-i/')s.T r ('+^ 

where in the last equality we have applied (pS)) . Let 

ri = max 775 t (28) 

0#SC/i0/TC/2 ' 

From (1271) and (12H1) we have that 

P(V,e/, ^.(^.) > 1 - 1/^, V,e/2 ^^(^^) > 1 - lA) 

(29) 

= <-^/''i.(i)E0^5C/,E0^TC/2(-l)"" + "^'^^^^T(O 

where As,T{t) = ^"(^/''■^'■^"^/'''l;^ ^(t) and L*^^ ^{t) = L^s .^(t)/L^{t) is a slowly varying function. 
Observe that, if S" C 5 and T' C T, then +00 > Ag' ;t'^) > As.rit) and, by the definition of 77, 
we have ^^^^(t) = 1 or Asj:-(t) as t — > 00, for all S C Ii and T C l2- Therefore, 

P{M{h) > 1 - l/t,M{h) > 1 - l/t) ^ t-^/'-iL^it). 

Moreover, considering rj = r/So,To for some Sq C /i,To C I2, and so Asg^Toit) = 1 < 
yieSo,j€ To, we must have A[iy^[jy{t) = 1, V i G 6*0,. 7 G Tq. Then 77 = 7?{i}.{j}, V i G ^o, j G Tq 
and ?7 < max^g/j j£72 But, by (|28p . 77 > maxig/j.jg/2 which leads to the result. □ 

In the following we present some examples where tail independence takes place. 
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Example 2.3. Consider {Ki}n>i o-n i-i-d. sequence of r.v.'s with distribution U(0, 1) and X = 
{Xi, X2, X3, X4) a random vector such that, Xi ~ min(V3, V2, Vi), X2 = niin(V4, V2, 14), ^3 = 
min(F4, V3,yi) and X4 ^ V5. Observe that, for < a; < 1, Fx^{x) = 1 - (1 - = Fx^ix) = 
Fx^ix) and Fx^ix) = x and hence F^^(x) = 1 ^ (1 ~ xY^'^ = F^^(x) — F^^(x) and Fj^^(x) = x. 
Consider Ii = {1,2} and I2 = {3,4}. 
We have successively, 

P{Fi{Xi) >l-t-\F3{X3) >l-t-^)=P{F2{X2) >l-t~\F3{X3) > =i-4/3^ 
and 

>l-t-\Fi{Xi) >l-t-^) = P{F2{X2) >l-t-\Fi{Xi) > =t-2. 
Hence, by Proposition \2.4\ we must derive f?({i.2},{3,4}) = 3/4. 
In fact, applying i27\). after some calculations we have 
P{M{h) > 1 - t-'x, Mih) > 1 - t-^y) 

j 2t-4/32.yi/3 ^ 2t-^xy - t-4/3a;4/3 - 2t-y^xy^/^ ~ ^t-^l^x^l'-y ,x<y 
\ t-*/3y2;i/3 + 2t-'^xy - 3t-'^/^x^/^y^ - t-'^/^x'^/^y , x > y. 

According to {24% coefficient rji^j-^ j^-^ can be obtained by taking x = y = 1 in the expression above, 
and by i23\) we obtain 



C({l,2},{3,4})(2;,y) 



22.yl/3 _ ^4/3 ^ X <y 

yx^l^ , X > y. 



which is homogeneous of order 4/3. 

Similarly, if we consider Ii = {1,2,3} and I2 = {4} we obtain ?7({i,2.3},{4}) = 1/2 ana 
C({i,2,3},{4})(a;, y) ^xy, and if Ii = {1} and I2 = {2,3,4} we have ?7({i},{2,3,4}) =3/4 and 

{xy^^^ , X < y 

2y^l/3_y4/3 =Cai,2},{3,4})(y,^). 



Example 2.4. Consider X = (^1, ■■■Xd) a standard d-variate Gaussian random vector with posi- 
tive definite correlation matrix. The bivariate tail- dependence structure is given by 

P{F,iX,) > 1 - l/t,Fj{Xj) > 1 - l/t) - Cp^^.t-2/(i+P.,.)(iog(t))-P.../(i+ft..), as t^ 00, (30) 

for i, j G {1, ...,d}, i < j , where pij = corr{Xi, Xj) ^ {—1,1} and 

Cp,,, = (1 + P,,jf^'{l - p.,,)"'/'(47r)-'''-/(i+'''-). 

Hence H25\} holds for I = {i} and J = {j} with ~ (l + jOij)/2 (see Ledford and Tawn J 12^ . 1996; 
Draisma et al. 2004)- According to Hua and Joe ffS^. 2011), i25fl also holds for non-empty 
sets /i,/2 C {l,...,d}. If we consider P(/i,/2) = inax{pij- : i £ Ii,j £ I2} then, by Proposition \2.4\ 
we find 'q(i-^j^) = (1 + /9(/j_/2))/2, provided the left-hand side of \24^ is non-null. 
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3 Estimation 

Several estimators for the bivariate stable tail dependence function in ([T7|) or even for the more 
general d-variate stable tail dependence function 

lim tp(Fi(Xi) > 1 - ^ V ... V F,i{Xd) >l-^) (31) 

have been considered in literature. For a survey, see Krajina (2010) [TT]. According to relation 
(fT5)) . they can be applied to our function /(^^'^^^(a;"^, y^^). 

We remark that these are based on asymptotic results that depend on a sequence of positive 
integers, {fc„}, going to infinity at a lower rate than n. For instance, the estimator based on p7p 
by plugging-in the respective empirical counterparts given by 

k fc 1 " 

n n 

where Fi{u) ~ n^^ X]fe=i '^{x^<u} is the empirical d.f. of / = 1, 2, is consistent and asymptoti- 
cally normal if {fc„} is an intermediate sequence, i.e., fc„ — > cxj and fc„/n — > 0, as n — >■ cx) (Huang 
1992 [S]). The choose of the value k in the sequence {fc„} that allows the better trade-off between 
bias and variance is of major difhculty, since small values of k come along with a large variance 
whenever an increasing k results in a strong bias. Therefore, simulation studies have been carried 
out in order to find the best value of k that allows this compromise. 

As mentioned before, the upper-tail dependence function in (jl4p can be viewed as an extension 
of the bivariate upper-tail dependence function of Schmidt and Stadtmiiller (2006) given in (fTO)) , 
by taking in this limit the random pair (A/(/i), M(/2)) instead of (F(Xi), F(X2)). The estima- 
tors considered in Schmidt and Stadtmiiller (2006), for which strong consistency and asymptotic 
normality have been established, allow to estimate our function A[|^'^^''(a;, y), as well as coefficient 
^i^j^) = A\j''^^\l, 1). However they are also based on asymptotic results with the same drawback 
of including an intermediate sequence, already referred above. 

In order to overcome this problem, we shall present a totally different and very simple approach. 
More precisely, the following result suggests an estimation procedure for the c?-variate stable tail 
dependence function in (j3ip that only evolves a sample mean. 

Proposition 3.1. Under the conditions of Provosition lKH we have, forl(xi, Xd) — — \ogF{xi, 

. ^ E{Fi{XiY^ y ...y FdjXdY^) 

""'^ 1 - E{F,{X^Y^ V ... V Fa{XaY-) ■ 
Dem. Consider for G{x) = exp(— l/x). Observe that 

E{G{xiXi) V ... V G{xdXa)) = E{G{xiXi V ... V xaXa)) 
and the d.f. of xiXi V ... V XdXd is given by 

P{xiXi V ... V XdXd < u) = P{Xi < u/xi, Xd < u/xd) 

= F{ux^\...,ux-') (32) 
= exp — l(uxi^ , ...,ux^^)) . 
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Hence 



cxp{—u^^) exp ( — l(^uxi^, ux^^)) — ( — l{^uxi^, mx^^)) 







_ _ _ _ d _ _ _ 

= J cxp{-u~^)exp{-u~'^l{x^'^,...,x2'^)) — {-u~^l{x^^,...,x';^^)) 



(33) 



l + ;(xr\...,a;^')' 

Now just observe that G{xiXi V ... V x^Xa) ^ FiiXi)^^""^ V ... V ^^(Xd)^/^''. □ 

Remark 3.1. Observe that the d-variate stable tail dependence function in iSl]) corresponds to 
-logFix^\...,x-'). 

By applying Proposition l3.1l with Xj replaced by xj^ , j = 1, d, we get the following corollary. 
Corollary 3.2. Under the conditions of Provosition we have 

^^^1-' >- i-E{M{hY/-y ^-^^^ 

yej,-l {y ,y > - i ^ EiM{hy/y) ^"^^^ 



and 



^ ' i-E{M{hY/^y M{hY/yy 



(36) 



Consider the estimators derived from Proposition 12.11 and Corollary 13.21 by plugging-in the 
respective sample means, respectively, 



^ ^ l-F,{X^r^y ...yPaiXaY^' ^ ' 



and 



= Jm!=, ye,^ = -^m^ and V^^^^^ i^-\y^^) = ^^(^^^^^ V ^^(^^)^^- . (33) 
l-M(/i)V^ l-M(/2)i/a l-A/(/i)V^V7\/(/2)i/y 

where 

^^^^^ = i E V ^.(^f )'/^ ^^^^c^ = ^ E V ^.(^i^V/^ (39) 

1=1 je/i 1=1 
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and 



1 " 

Mihy/-\/ Mihy/y = -J2i\/ Fjixff'- V V F,{X. 



n 

1=1 jell jeh 



(40) 



We will consider two situations: the first one for known margins and the second one for unknown 
margins. 

In case the margins are known, they become unit Frechet by transformation —1/ log Fj(Xj) for 

J el c{i,...,d}. 

It is quite straightforward to deduce the consistency and asymptotic normality of estimators 
([57)1 and dnni) by the well-known Delta Method. 

Proposition 3.3. Under the conditions of Provosition WA[ we have 

V^(l{xi,...,Xd)-l{xi,...,Xd)) ^ N{0,a'^), (41) 

where l{xi, Xd) is the estimator derived from Provosition \3.1\ by plugging-in the respective sample 
mean given in i37\) and 

2 _ l{xi,...,Xd}(^l+l{xi,...,Xd)) 

{2+l{xi,...,xa)) ' 

Dem. Let Yi, i = 1, n, be independent copies of F = Fx{X^^^ \l Fd{Xd)^''- We have 
that V^(F-^y) ^ N{0,a^.), where hy = E{Fi{XiY^\J ...\J FdiXdY"") and = Var{Fi{XiY^y 
... V Fd^XdY"^)- By a similar reasoning of (|33|) we derive 

E{{F,{X,Y^ V ... V Fd[XdY^f) - sfe^ 

and hence, 

Var{{F,{XiT' V ... V Fd{XdY'?) = 7 ^1-1^ 

[2+l(xx,...,xa)) [l+l{xi,...,xa)) 

Let g{x) = (1 - xY^ - 1- We have [^'(/iy)]^ = (1 - ^jy)"** and, by the Delta Method, s/n{g(y) - 
xei,)^ N{Q, alii- tivY^)- □ 

Corollary 3.4. Under the conditions of Provosition lKTX we have 

V^ixei, - xei, ) -> iV(0, aj), (42) 

V^iyei, - yei^) -> iV(0, al) (43) 

and 

V^{T^'^^'^\x-\y-') -l^'^-'^\x'\y-')) ^ N{0,al), (44) 
where xe/^ , y'ej^ and l^^^'^^\x~^ ^y~^) are given in VS^) and 



2 X€l^{l + XtlY) 
(Ti — 



(2 + xe/J ' 



2 _ ye/2(i + ye/2)^ 



(2 + yti,) 

V 2 



{2 + lii^^i--){x-\y-^)) 



(45) 
(46) 

(47) 
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Based on the definition in (|14l) . a natural estimator for the upper-tail dependence function is 

A[|-^^)(x,2/) =xe7: +2/?}; ^V^^'^\x-\y-'), (48) 

with xeYn ^iid l'^^^'^^\x~^ , y~^) stated in ([55|) . Hence we have the following estimator for the 
extremal coefficient of dependence between X/^ and X/^: 

Hh.h) = + e}; - ?7iu/2 ■ (49) 

where ?/,u/2 = I^^^^^^' (1, !)■ 

Proposition 3.5. Estimators l{xi, ...,Xd) and A[j^''^^\x,y) in i37\) and i48\ ), respectively, are 
strong consistent. Consequently, the same holds for in ^9\ l. 

Dem. Just observe that, as the sample mean Af(/i)i/^ converges almost surely to the 
mean value £;(M(Ji)i/^), i.e., M{hy/^ ^ E{M{h)^^''), then^xeT^ = g{JlihW^)J^ xei, = 
g{E{M{Iiy/-^)), where g{x) = (l-x)"i-l. Analogously for ye^, rt^i'-f^)^^,-!^ ^-i) and 7(2:1, Xd). 
Now, the strong consistency of A^^''^^\x, y) is straightforward from 

\V-^''''\x,y) - K^i''''\x,y)\ < \x^, - xej,\ + \y^, - yej,\ + \li'^-i-\x-\y-') - l('^''-) {x-\y-% □ 



Now consider Fj the empirical d.f. of Fj, j ~ l,...,d, 

1 " 



k=l 



The denominator n + 1 instead of n in the empirical d.f. concerns estimation accuracy and other 
modifications can be used. For a discussion see, for instance, Beirlant et al. [T] (2004). 

In case of unknown margins, we can replace Fj by the respective empirical d.f. Fj, j = 1, d, 
in ([57)1 and ([55]) . More precisely, we have 



l{xi,...,Xd) = _ =, (50) 



1 - FiiXiY^ W ... W FdiXdY^ 

as well as, 



xt,^ = _JIM=, . and P^^^^\x-\y-^) ^ ^^hY''^ ^ m.)^'^ (51) 



where 



— 1 " 

Fi(Xi)-i V ... V Fd{Xd)-^ V 



n 

i=l jG{l,...,d} 



1=1 je/i 1=1 jG/2 



(52) 



(53) 



and 



, n 

M{hY/^ V MihY/y - - E ( V F^i^i'^y^" V V F,(xf )i/^). (54) 



n 



We still have asymptotic normality of estimators in from the following result stated 

in Fermanian et al. (2002, [H]. Theorem 6). 
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Theorem 3.6. (Fermanian et al. (2002) Theorem 6) Let F have continuous marginals and let 
copula Cp in |0i have continuous partial derivatives. Then 



1 " r 
^^{J(i?i(xW),...,P,(xW))_i?(J(Fi(X«),...,F,(X«)))}^ / 



G(wi, Ud)dJ{ui, Ud) 

[04]'' 



in distribution in ^°°([0, 1]''), where the limiting process and G are centered Gaussian, and J : 
[0, 1]'' — ^ R is of hounded variation, continuous from above and with discontinuities of the first kind 
(Neuhaus, 1971 pf). 

The asymptotic normality of estimators (j50p and (jSip is now derived from a general version of 
the Delta Method as considered in Schmidt and Stadtmiiller [52] (2006; Theorem 13). 

We also state strong consistency of estimators l{xi, .■.,Xd) in ([50]) and 

Vi''''\x,y)^xrr,+yrr.-^''-''\x-\y-'), (55) 

with Xff/j, ye]~2 and /^^^'^^^(a;^^, y^^) given in ([5T|) . and hence of estimator 

^(/i J2) = + - e/iU72 7 (56) 

where e/iu/2 = i^^^-^"' (1. !)■ 

Proposition 3.7. Estimators l{xi, ...,Xd) in 15(J\) and A^'''^^\x,y) in !i55\) are strong consistent. 
Therefore, the same holds for estimator ^j^ j,^^ in i56]} . 

Dem. The proof runs along the same lines as the one of Proposition 13.51 We only prove the 
more general case l(xi, Xd) —4 l{xi, Xd)- Observe that 

^ I^E V ^(^fr-^E V ^.(^fr^- 

i=i ]e{i,...,d} 1=1 je{i,...,d} 

+|^E V F,{xfr-E{ V F^ix^r), 

■'=iie{i d} je{i d} 

where the second term converges almost surely to zero by the Strong Law of Large Numbers. 
For the first term we have, successively, 

1^ ElLi Vje{i,...,<i} FjiX^j'^y-^^ - i J2'i=i Vje{i,...,rf} FjiX^'^)"^^ I 

<^E:uv,e{i d}\FAx!f^r'-F,{xfr^\ 

which converges almost surely to zero according to Gilat and Hill ([7], 1992; proof of Theorem 1.1). 
□ 
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4 Application to financial data 

In this section we show that tail dependence is present in financial data. Our analysis is based 
on negative log-returns of daily closing values of the stock market indexes, CAC 40 (France), 
FTSEIOO (UK), SMI (Swiss), XDAX (German), Dow Jones (USA), Nasdaq (USA), SP500 (USA), 
HSI (China), Nikkei (Japan). The period covered is January 1993 to March 2004. More precisely, 
we consider the monthly maximums in each market and group the indexes in Europe (CAC 40, 
FTSEIOO, SMI, XDAX), USA (Dow Jones, Nasdaq) and Far East (HSI, Nikkei). The scatter plots 
in Figure [T] show the presence of dependence between the monthly maximums in Europe and USA, 
Europe and Far East, USA and Far East, respectively. We are interested in assessing the amount 



— I 1 1 1 r 

0.2 0.4 0.6 0.8 1.0 



1 1 1 1 

0.2 0.4 0.6 0.8 1.0 




0.2 0.4 0.6 0.8 
USA 



Figure 1: Scatter plots of the monthly maximums (84 data points) in Europe versus USA, Europe 
versus Far East and USA versus Far East. 



of tail dependence between the three big world markets referred: Europe, USA and Far East, and 
this can be achieved through the extremal coefficient of dependence ^(i-ij^)^ defined in (|15p . As we 
do not know the margins distribution, we use estimator e(/j j,) in (j56p based on ranks. In Table 
[l]are the obtained estimates for several groups, Ii and l2- One can see that the Far East market 
has less infiuence (lower values of the coefficient) but Europe and USA have a stronger effect on 
each other and on the respective group of foreign markets. Observe that the difference between 
these two magnitudes of dependence is almost in the proportion 1:2. 



h 


h 




Europe 


USA 


1.008324625 


Europe 


Far East 


0.568780467 


USA 


Far East 


0.364358832 


Europe 


USA U Far East 


1.125919957 


USA 


Europe U Far East 


0.921498322 


Far East 


USA U Europe 


0.481954164 



Table 1: Estimates of the extremal coefficient of dependence e^/ijj) ^'^^ ^^-^^ indicated groups, Ii 
and 12- 



5 Conclusion 

In this work we introduce a new upper-tail dependence concept for a random vector which extends 
the one in Schmidt and Stadtmiiller ([52], 2006). Our approach weakens the usual imposed multi- 
variate tail dependence and can be treated with bivariate techniques. The new function extends 
the well-known relation of Huang ([2] 1992) for a MEV with unit Frechet marginals, and gives rise 
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to the so-called extremal coejficient of dependence as it is expressed through the extremal coeffi- 
cient in Tiago de Oliveira ([2B], 1962-63) and Smith ([131j 1990). We also enlarge our discussion 
to tail independence in the sense of Ledford and Tawn ([12] [121 j 1996, 1997). At this point we are 
beyond MEV distributions which only admit tail dependence or (exact) independence. 

In calculating the moments of the r.v.'s involved in our function, we arrive at very simple 
estimators whose asymptotic normality is stated. These can also be applied to the well-known 
stable tail dependence function. We also prove strong consistency of the proposed estimators for 
our measures. We end with an application to financial data presenting tail dependence. 
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